home *** CD-ROM | disk | FTP | other *** search
/ BCI NET 2 / BCI NET 2.iso / archives / misc / integra1_0a.lha / Integra / examples / 1test.c < prev    next >
Encoding:
C/C++ Source or Header  |  1994-10-22  |  1.1 KB  |  51 lines
  1. /*
  2.     Test per Int_eq_points
  3.  
  4. */
  5. #include "simple.h"
  6. #include <stdio.h>
  7. #include <stdlib.h>
  8. #include <math.h>
  9.  
  10.  
  11. int main(int argc, char *argv[])
  12. {
  13. int i, n;
  14. double d, x, fx[1000];
  15.  
  16. if (argc==1) {
  17.     printf("Quanti punti? ");
  18.     scanf("%d",&n);
  19. } else n=atoi(argv[1]);
  20. if (n>1000) n=1000;
  21. if (n<2) n=2;
  22.  
  23. printf("Test di Int_eq_points per la funzione x tra 2 e 5\n");
  24. for (i=0; i<n; i++) {
  25.     x = 2 + i*(5-2)/(double)(n-1);
  26.     fx[i] = x;         /*  x*x*x*(x-1); */
  27. }
  28. d = Int_eq_points(n, (5-2)/(double)(n-1), fx);
  29. printf("Esatto    : 10.5\n");
  30. printf("%-2d punti  : %.16g\n",n,d);
  31.  
  32. printf("Test di Int_eq_points per la funzione x^4-x^3 tra 2 e 5\n");
  33. for (i=0; i<n; i++) {
  34.     x = 2 + i*(5-2)/(double)(n-1);
  35.     fx[i] = x*x*x*(x-1);
  36. }
  37. d = Int_eq_points(n, (5-2)/(double)(n-1), fx);
  38. printf("Esatto    : 466.35\n");
  39. printf("%-2d punti  : %.16g\n",n,d);
  40.  
  41. printf("Test di Int_eq_points per la funzione x*sin(x)^2 tra 2 e 5\n");
  42. for (i=0; i<n; i++) {
  43.     x = 2 + i*(5-2)/(double)(n-1);
  44.     fx[i] = x*sin(x)*sin(x);
  45. }
  46. d = Int_eq_points(n, (5-2)/(double)(n-1), fx);
  47. printf("Esatto    : 5.5748036294843\n");
  48. printf("%-2d punti  : %.16g\n",n,d);
  49.  
  50. }
  51.